home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_nrm2.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  2.8 KB  |  144 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_nrm2 (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int N = 1;
  13.    float X[] = { 0.317f };
  14.    int incX = -1;
  15.    float expected = 0.0f;
  16.    float f;
  17.    f = cblas_snrm2(N, X, incX);
  18.    gsl_test_rel(f, expected, flteps, "snrm2(case 28)");
  19.   };
  20.  
  21.  
  22.   {
  23.    int N = 1;
  24.    double X[] = { 0.071 };
  25.    int incX = -1;
  26.    double expected = 0;
  27.    double f;
  28.    f = cblas_dnrm2(N, X, incX);
  29.    gsl_test_rel(f, expected, dbleps, "dnrm2(case 29)");
  30.   };
  31.  
  32.  
  33.   {
  34.    int N = 1;
  35.    float X[] = { 0.776f, 0.983f };
  36.    int incX = -1;
  37.    float expected = 0.0f;
  38.    float f;
  39.    f = cblas_scnrm2(N, X, incX);
  40.    gsl_test_rel(f, expected, flteps, "scnrm2(case 30)");
  41.   };
  42.  
  43.  
  44.   {
  45.    int N = 1;
  46.    double X[] = { 0.549, -0.354 };
  47.    int incX = -1;
  48.    double expected = 0;
  49.    double f;
  50.    f = cblas_dznrm2(N, X, incX);
  51.    gsl_test_rel(f, expected, dbleps, "dznrm2(case 31)");
  52.   };
  53.  
  54.  
  55.   {
  56.    int N = 2;
  57.    float X[] = { 0.14f, -0.632f };
  58.    int incX = 1;
  59.    float expected = 0.647320631527f;
  60.    float f;
  61.    f = cblas_snrm2(N, X, incX);
  62.    gsl_test_rel(f, expected, flteps, "snrm2(case 32)");
  63.   };
  64.  
  65.  
  66.   {
  67.    int N = 2;
  68.    double X[] = { 0.696, -0.804 };
  69.    int incX = 1;
  70.    double expected = 1.06340584915;
  71.    double f;
  72.    f = cblas_dnrm2(N, X, incX);
  73.    gsl_test_rel(f, expected, dbleps, "dnrm2(case 33)");
  74.   };
  75.  
  76.  
  77.   {
  78.    int N = 2;
  79.    float X[] = { 0.281f, -0.063f, 0.367f, 0.232f };
  80.    int incX = 1;
  81.    float expected = 0.521001919382f;
  82.    float f;
  83.    f = cblas_scnrm2(N, X, incX);
  84.    gsl_test_rel(f, expected, flteps, "scnrm2(case 34)");
  85.   };
  86.  
  87.  
  88.   {
  89.    int N = 2;
  90.    double X[] = { -0.359, -0.76, -0.906, -0.108 };
  91.    int incX = 1;
  92.    double expected = 1.24055672986;
  93.    double f;
  94.    f = cblas_dznrm2(N, X, incX);
  95.    gsl_test_rel(f, expected, dbleps, "dznrm2(case 35)");
  96.   };
  97.  
  98.  
  99.   {
  100.    int N = 2;
  101.    float X[] = { 0.918f, -0.126f };
  102.    int incX = -1;
  103.    float expected = 0.0f;
  104.    float f;
  105.    f = cblas_snrm2(N, X, incX);
  106.    gsl_test_rel(f, expected, flteps, "snrm2(case 36)");
  107.   };
  108.  
  109.  
  110.   {
  111.    int N = 2;
  112.    double X[] = { 0.217, -0.588 };
  113.    int incX = -1;
  114.    double expected = 0;
  115.    double f;
  116.    f = cblas_dnrm2(N, X, incX);
  117.    gsl_test_rel(f, expected, dbleps, "dnrm2(case 37)");
  118.   };
  119.  
  120.  
  121.   {
  122.    int N = 2;
  123.    float X[] = { 0.31f, 0.059f, -0.442f, 0.987f };
  124.    int incX = -1;
  125.    float expected = 0.0f;
  126.    float f;
  127.    f = cblas_scnrm2(N, X, incX);
  128.    gsl_test_rel(f, expected, flteps, "scnrm2(case 38)");
  129.   };
  130.  
  131.  
  132.   {
  133.    int N = 2;
  134.    double X[] = { 0.609, 0.615, -0.143, -0.957 };
  135.    int incX = -1;
  136.    double expected = 0;
  137.    double f;
  138.    f = cblas_dznrm2(N, X, incX);
  139.    gsl_test_rel(f, expected, dbleps, "dznrm2(case 39)");
  140.   };
  141.  
  142.  
  143. }
  144.